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Abstract 

Background: A bacterial strain previously isolated from pyrite mine drainage and named BAS-10 was tentatively 
identified as Klebsiella oxytoca. Unlikely other enterobacteria, BAS-10 is able to grow on Fe(lll)-citrate as sole carbon 
and energy source, yielding acetic acid and C0 2 coupled with Fe(lll) reduction to Fe(ll) and showing unusual 
physiological characteristics. In fact, under this growth condition, BAS-10 produces an exopolysaccharide (EPS) 
having a high rhamnose content and metal-binding properties, whose biotechnological applications were proven 
as very relevant. 

Results: Further phylogenetic analysis, based on 16S rDNA sequence, definitively confirmed that BAS-10 belongs to 
K. oxytoca species. In order to rationalize the biochemical peculiarities of this unusual enterobacteriun, combined 
2D-Differential Gel Electrophoresis (2D-DIGE) analysis and mass spectrometry procedures were used to investigate 
its proteomic changes: i) under aerobic or anaerobic cultivation with Fe(lll)-citrate as sole carbon source; ii) under 
anaerobic cultivations using Na(l)-citrate or Fe(lll)-citrate as sole carbon source. Combining data from these 
differential studies peculiar levels of outer membrane proteins, key regulatory factors of carbon and nitrogen 
metabolism and enzymes involved in TCA cycle and sugar biosynthesis or required for citrate fermentation and 
stress response during anaerobic growth on Fe(lll)-citrate were revealed. The protein differential regulation seems to 
ensure efficient cell growth coupled with EPS production by adapting metabolic and biochemical processes in 
order to face iron toxicity and to optimize energy production. 

Conclusion: Differential proteomics provided insights on the molecular mechanisms necessary for anaeorobic 
utilization of Fe(lll)-citrate in a biotechnologically promising enterobacteriun, also revealing genes that can be 
targeted for the rational design of high-yielding EPS producer strains. 




MICROBIAL CELL 
FACTORIES 



Background 

Several species of enterobacteria use citrate as sole car- 
bon and energy source. This capability is firstly due to 
appropriate transporters for citrate up-take, such as the 
citrate-specific proteins CitH and CitS [1,2] or like the 
tripartite tricarboxylate transporter (TTT) TctABC sys- 
tem able to transport several tricarboxylic acids into the 
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bacterial cell [3] or like the ferric citrate transport sys- 
tem (the product of the fecABCDE operon) that shuttles 
the Fe (III) -citrate complex into the cytoplasm [4-7]. 

During aerobiosis, intracellular citrate is catabolized 
throughout the TCA cycle. Under anaerobic conditions, 
when TCA cycle is down- regulated, enterobacteria spe- 
cies, like Klebsiella pneumoniae and Salmonella typhi- 
murium, can grow on citrate by a Na(I) -dependent 
pathway, forming acetic acid and C0 2 as final metabo- 
lites [1,8,9]. Genes specific for anaerobic citrate fermen- 
tation, such as those coding for regulators (citAB), 
catabolic enzymes (citCDEFG and oadGAB) and citrate 
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transporters (citS and citW), have been identified in 
these bacteria [1,2,8,10]. The presence of sodium is es- 
sential for citrate symport by CitS and for the activity of 
oadGAB gene products that form the oxaloacetate de- 
carboxylase complex. In fact, oxaloacetate decarboxylase 
converts oxaloacetate into pyruvate and pumps sodium 
externally to synthesize ATP [1,2,8,9]. 
Generally, iron is one of the major limiting nutrients 

[II] and citrate-fermenting enterobacteria do not usually 
thrive on high concentrations of Fe(III)-citrate as sole 
carbon and energy source [12]. Indeed, there are habitats 
where the abundance of Fe(III) is so high, like in pyrite 
mine drainages, which represents one of the major ele- 
ments to make rust-red the acidic waters. In this case, 
iron can represent an environmental hazard for life, es- 
pecially for its oxidative properties and for the presence 
of other metals which increase the total toxicity of mine 
drainages and cause a significant reduction of microbial 
biodiversity [13-15]. Only specialized species can survive 
in extreme habitats with high heavy metal concentra- 
tions and carbon-depleted conditions and Enter obactera- 
ceae are not expected to survive in such environments. 
Nevertheless, an enterobacterial strain was isolated 
under an iron mat formed by waters leached from pyr- 
ite mine drainages of Colline Metallifere, Tuscany, Italy 
[16]. This isolate, named BAS-10, was tentatively iden- 
tified as K. oxytoca on the basis of partial (422 nt) 16S 
rDNA sequence and API Enterotube test [16]. Unique 
among Klebsiella strains, BAS-10 can ferment and pro- 
liferate on Fe(III)-citrate as sole carbon and energy 
source, forming acetic acid and C0 2 coupled with Fe 

(III) reduction to Fe(II) [12]. Under these growth con- 
ditions, BAS-10 produces an EPS made of rhamnose 
(57.1%), glucuronic acid (28.6%) and galactose (14.3%), 
which shows metal-binding properties [17,18]. Al- 
though extracellular EPS have been reported over re- 
cent decades and their composition, structure, 
biosynthesis and functional properties have been exten- 
sively studied, only a few have been industrially devel- 
oped [19]. 

Chelating and sugar compositional properties of 
BAS-10 EPS are of high interest for potential biomed- 
ical, food, and environmental applications [18-20]. Fur- 
ther studies on BAS-10 physiology may be useful to 
develop efficient fermentation processes for EPS pro- 
duction. In order to investigate iron-dependent bio- 
chemical and metabolic adaptations and regulatory 
networks thereof during anaerobic growth on Fe(III)- 
citrate as the sole carbon source, a differential prote- 
omic approach was used. In particular, proteomic 
repertoires from BAS-10 grown on Fe(III)-citrate under 
anaerobiosis, on Na-citrate under anaerobiosis and on 
Fe (III) -citrate under aerobiosis were comparatively 
evaluated. 



Results and discussion 

Phylogenetic identification 

Physiological studies [12], EPS composition and metal- 
binding activity thereof [17,18] revealed characteristic 
peculiarities of BAS-10 strain. Thus, a sequence of 1447 nt 
gene was generated from BAS-10 16S rDNA (Add- 
itional File 1) to perform phylogenetic clustering. Two 
ClustalW analyses were performed by using the first 
twenty hits from Blast analysis, selecting only cultivable 
reference or whole strains from NCBI database (http:// 
blast.ncbi.nlm.nih.gov/) [21], respectively. Phylogenetic 
threes revealed that BAS-10 is related to K. oxytoca 
reference strain (Figure 1) and can be clustered into a 
group comprising other five very close related strains 
classified as K. oxytoca (Additional file 1 Figure SI). 
This investigation definitively confirmed previous taxo- 
nomical observations on BAS-10 [16]. 

Proteomic analysis 

To investigate iron-dependent cell processes and regula- 
tory networks thereof associated to citrate fermentation 
and EPS production, BAS-10 was cultivated: i) under 
aerobic conditions on Fe(III)-citrate (FEC) as sole car- 
bon source; ii) under anaerobic conditions using FEC or 
iii) Na(I)-citrate (NAC) as sole carbon source. Then, 
bacterial proteomes were analyzed using anaerobic FEC 
growth as pivotal condition. In particular, the "anaerobic 
FEC Vs aerobic FEC" differential analysis was carried 
out to reveal proteins that are required for anaerobic 
growth on FEC while "anaerobic FEC Vs anaerobic 
NAC" differential analysis was carried out to reveal pro- 
teins whose abundance is iron-dependent during anaer- 
obic growth. In all cases, BAS-10 biomass samples were 
all harvested at late exponential growth phases, when 
anaerobic growth on FEC is coupled with iron-binding 
EPS production [12,16] at an average yield of 6.7 g/L 
(dry weight EPS per culture volume). Proteomic 2D- 
DIGE maps showed the occurrence of 869 reproducible 
protein spots (i. e. detected in all 2D-Gels) which were 
used to create the match set for comparative analysis. 
According to the criteria reported in the experimental 
section for identification of deregulated proteins, com- 
parative analysis between FEC medium cultivations 
under anaerobic or aerobic conditions revealed 27 and 
32 protein spots up- and down-regulated under anaer- 
obic conditions, respectively (Table 1; Additional file 1 
Figure S2). On the other hand, comparative analysis be- 
tween cultivations in FEC and NAC media under anaer- 
obic conditions revealed 28 and 42 protein spots as up- 
and down-regulated in the presence of Fe(III) (Table 1; 
Additional file 1 Figure S2). Among these protein spots, 
46 resulted differentially abundant with a concordant 
profile in both analyses (anaerobic FEC vs aerobic FEC 
and anaerobic FEC vs anaerobic NAC) (Table 1), thus 
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Raoultella terrigena 84 -Y17658 
Enterobacter amnigenus JCM1237 -AB004749 
Yokenella regensburgei JCM2403 -GQ856082 
Enterobacter kobei CIP1 05566 -AJ508301 
Klebsiella oxytoca ATCC13182 -AF1 29440 
Klebsiella oxytoca BAS-10 
Citrobacter gillenii CDC4693-86 -AF025367 
Enterobacter nimipressuralis LMG10245 -Z96077 
j Enterobacter ludwigii EN-1 19 -AJ853891 
V- Enterobacter cancerogenus LMG2693 -Z96078 
L Leclercia adecarboxylata LMG2803 -GQ856082 
- Enterobacter hormaechei ATCC49162 -AJ508302 
Escherichia vulneris ATCC33821 -AF530476 
- Enterobacter cloacae 279-56 -AJ251469 
— Raoultella ornithinolytica CIP1 03364 -U78182 
Streptomyces tendae ATCC19812 -D63873 



Figure 1 Phylogenetic analysis of BAS-10 strain performed by using 16S rDNA sequences. The first twenty hits from BLAST analysis 
performed by selecting only cultivable reference strains were used to generate the phylogenetic three. The 16S rDNA sequence of Streptomyces 
tendae was used as outgroup. Distance unit is based on sequence identity. NCBI accession number of each 16S rDNA is reported after hyphen. 



indicating that these proteins are specifically associated 
with anaerobic growth on FEC The selected differen- 
tially abundant protein spots (83 in number) were alky- 
lated, trypsinolyzed and identified by MALDI-TOF- 
peptide mass fingerprinting (PMF) and nanoLC-ESI- 
LIT-MS/MS procedures (Table 1 and Additional file 1 
Table SI). Corresponding protein entries were clustered 
into functional groups, according to KEGG (http://www. 
genome.ad.jp/kegg/kegg2.html) [22], EcoCyc (http:// 
ecocyc.org/) [23] and ExPASy (http://www.expasy.org) 
[24] databases (Additional file 1 Figure S3). Among 
the identified proteins, the most represented ones were 
related to central carbon metabolism -including gly- 
colysis/gluconeogenesis, tricarboxylic acid (TCA) cycle, 
citrate fermentation, pentose phosphate (PP) pathway 
-and membrane transport (Table 1; Additional file 1 
Figure S3). 

Proteins required for Fe(lll)-citrate fermentation 

The differential proteomic analysis revealed up- 
regulation of citrate lyase a- and (3-subunit (CitE and 
CitF), oxaloacetate decarboxylase a-subunit (OadA), 
pyruvate formate lyase (PFL), phosphotrans-acetylase 
(PAT) and acetate/propionate kinase (ACK) in anaerobic 
FEC with respect to both aerobic FEC and anaerobic 
NAC (Table 1; Figures 2 and 3). These proteins are 
known to be involved in anaerobic citrate fermentation, 
eventually converting citrate to acetate with production 
of ATP. In particular, CitF and CitE are part of the cit- 
rate lyase complex, the key enzyme in initiating the 



anaerobic utilization of citrate, which is responsible for 
catalyzing the conversion of citrate into acetate and oxa- 
loacetate. The oxaloacetate decarboxylase complex, con- 
stituted by alpha, beta and gamma subunits, catalyzes 
the second step of citrate fermentation by converting 
oxaloacetate into pyruvate, the latter being the substrate 
of PFL, a typical enzyme of enterobacteria growing 
under anaerobic conditions [25], that generates formate 
and acetyl-CoA. Phosphotrans-acetylase (PAT) catalyzes 
conversion of acetyl-CoA into acetyl-P, which finally 
transfers the phosphate group to ADP to yield acetic 
acid and ATP in the ACK-catalyzed last step of citrate 
fermentation. Different forms of CitE, CitF, OadA and 
PFL were identified in K. oxytoca BAS-10 2D-protein 
maps (Table 1, Figures 2 and 3; Additional file 1 Table 
SI and Figure S2). The occurrence of several protein 
forms having Mw and/or pi different from the predicted 
one would imply that protein activity could be con- 
trolled by post-translational modifications (PTMs), like 
covalently-bound charged small molecules (as in the 
case horizontal spot trains) or proteolytic digestion (as 
in the case of protein isoform with reduced Mw). In- 
deed, PTMs, like acetylation, succinylation and radical 
formation have already been shown for PFL activity 
regulation [26-28], thus giving count for the generation 
of different protein spots in 2D-maps. In the case of 
OadA, biotinylation in lysine, reported to be crucial for 
oxaloacetate decarboxylase activity [29], could explain 
the slight increase of the measured Mw while the gener- 
ation of reduced Mw species is likely to be due to 
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Table 1 Functional classification and relative fold change of K. oxytoca BAS-10 proteins differentially expressed during 
anaerobic growth in FEC 



Spot Protein name 



Acronym 



NCBI 
code 



Theor. 

pi/Mr 

(kDa) 



Exp. 

pi/Mr 

(kDa) 



Functional 
class 3 



Relative fold Relative fold 
change b (-0 2 change b (-0 2 
FEC Vs +0 2 FEC) FEC Vs -0 2 NAC) 



Glc 



1 NADP-specific glutamate GLDH 376386651 6.73/49 6.26/45 
dehydrogenase 

5* Nitrogen regulatory GlnK 376394829 5.84/12 5.62/13 
protein P-ll 2 

6 Nitrogen regulatory GlnK 376394829 5.84/12 5.62/13 
protein P-ll 2 

7 Nitrogen regulatory GlnK 376394829 5.84/12 5.20/13 
protein P-ll 2 

3* Ketol-acid reductoisomerase llvC 376396684 5.25/54 5.04/54 



4 D-3-phosphoglycerate PHGDH 365906793 6.06/44 5.92/46 
dehydrogenase 

8 Major outer membrane Lpp 376401994 9.36/8 4.16/12 
lipoprotein 1 

9 Major outer membrane Lpp 376401995 9.36/9 4.04/12 
lipoprotein 1 

10 Major outer membrane Lpp 376401994 9.36/8 4.44/12 
lipoprotein 1 

37 Aerobic respiration control ARCA 376394332 5.30/27 5.04/29 
protein ArcA 

19* Glucose-specific EIIA Glc 376399349 4.73/18 4.6/22 
phosphotransferase 
enzyme HA component 

20* Glucose-specific EIIA' 
phosphotransferase 
enzyme HA component 

17* Hypothetical protein SucB 
HMPREF9694_04828 
(dihyd rolipoyllysine-residue 
succinyltransferase, 
E2 component) 

41 Succinyl-CoA synthetase A-SCS 365909192 5.89/30 5.65/30 
subunit alpha 

33 Succinyl-CoA synthetase B-SCS 365909191 5.35/42 5.05/43 
subunit beta 

42 Fumarate hydratase class II FH 376401593 6.02/50 5.92/46 

2* Fumarate reductase FrdA 376396036 5.60/66 5.44/77 

flavoprotein subunit 

22* Malate dehydrogenase MDH 376398117 5.57/33 5.34/34 

23* Malate dehydrogenase MDH 376398117 5.57/33 5.08/29 

14 Citrate lyase beta subunit CitE 376376734 5.33/31 4.87/29 

26 Citrate lyase beta subunit CitE 376376734 5.33/31 5.07/31 

36* Citrate lyase beta subunit CitE 376376734 5.33/31 5.19/29 



376399349 4.73/18 4.62/20 



376395176 5.74/44 5.39/54 



Amino Acid 
Metabolism 

Amino Acid 
Metabolism 

Amino Acid 
Metabolism 

Amino Acid 
Metabolism 

Amino Acid 
Metabolism 

Amino Acid 
Metabolism 

Cell wall 



Cell wall 

Cell wall 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 



Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 

Central carbon 
metabolism 



1 

-3.2 

1 

1 

2.3 

1 

1 

1 

1 

1.9 
-1.9 

-1.5 



-3.9 
-1.8 
1 

2.1 

-3.1 

-1.9 

1.9 

4.9 

3.0 



2.6 

-8.7 

-2.7 

-2.2 

2.2 

-2.7 

3.4 

3.4 

1.7 

1 

-1.7 
-1.4 
-3.6 

1 
1 

-2.7 

2.6 

-2.7 

-2.7 

1 

1 

2.4 
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Table 1 Functional classification and relative fold change of K. oxytoca BAS-10 proteins differentially expressed during 
anaerobic growth in FEC (Continued) 



15 Citrate lyase, alpha subunit CitF 365908401 5.94/55 

16* Citrate lyase, alpha subunit CitF 365908401 5.94/55 

24* Oxaloacetate decarboxylase OadA 376376739 5.44/63 
alpha chain 

25* Oxaloacetate decarboxylase OadA 376376739 5.44/63 
alpha chain 

81* Oxaloacetate decarboxylase OadA 376376739 5.44/63 
alpha chain 

29* Pyruvate formate lyase PFL 376393708 5.63/85 

30* Pyruvate formate lyase PFL 376393708 5.63/85 

31* Pyruvate formate lyase PFL 376393708 5.63/85 

32* Pyruvate formate lyase PFL 376393708 5.63/85 

35* Pyruvate formate lyase PFL 376393708 5.63/85 

27* Phosphotrans-acetylase PTA 365911544 5.26/77 

1 2* Acetate kinase A/propionate ACK 36591 1 543 5.84/44 

kinase 2 

46* Triosephosphate isomerase TIM 376397984 5.82/26 

38* Pyruvate kinase PK 365911037 6.00/52 

21* Glyceraldehyde-3-phosphate GAPDH 376401359 6.33/36 
dehydrogenase 

28 Pyruvate dehydrogenase PDH 365908482 5.47/100 

subunit E1 

40 3-Methyl-2-oxobutanoate PanB 365908516 5.64/28 
hydroxymethyltransferase 

49 6,7-Dimethyl-8-ribityllumazine RibH 365908789 5.12/16 
synthase 

51 DNA protection during DPS 376393609 5.72/19 

starvation protein 

50 DNA-binding protein HU-alpha HU-2 376396790 9.40/9 

52* ATP synthase subunit beta ATP(3 376397178 4.93/50 

53* ABC transporter arginine- ArtJ 376393664 6.90/27 

binding protein 1 

57* Glutamate and aspartate DEBP 365909131 8.61/33 
transporter subunit 

58* Glutamine-binding periplasmic GlnH 376388520 8.74/27 
protein 

59* Glutamine-binding periplasmic GlnH 376388520 8.74/27 
protein 

60* Glutamine-binding periplasmic GlnH 376388520 8.74/27 
protein 



5.98/55 Central carbon 8.7 
metabolism 

5.89/55 Central carbon 14.1 
metabolism 

5.04/64 Central carbon 10.4 
metabolism 

4.7/43 Central carbon 3.7 
metabolism 

4.74/43 Central carbon 4.0 
metabolism 

5.09/80 Central carbon 2.7 
metabolism 

5.3/80 Central carbon 2.1 
metabolism 

5.37/80 Central carbon 2.9 
metabolism 

5.44/77 Central carbon 2.7 
metabolism 

5.16/43 Central carbon 2.3 
metabolism 

5.27/65 Central carbon 3.6 
metabolism 

5.39/27 Central carbon 4.1 
metabolism 

5.5/26 Central carbon 1.6 
metabolism 

5.78/55 Central carbon 1.9 
metabolism 

6.2/31 Central carbon 2.2 
metabolism 

5.52/67 Central carbon 2.4 
metabolism 

5.34/29 Cofactor biosynthesis 1 

4.87/14 Cofactor biosynthesis 1 

5.61/16 DNA ricombination, -1.5 
replication and repair 

9.07/1 1 DNA ricombination, 1 
replication and repair 

4.81/51 Energy metabolism 2.5 

5.98/25 Membrane Transport 4.9 

7.85/30 Membrane Transport -7.7 

5.81/26 Membrane Transport -1.9 

7.8/25 Membrane Transport -3.1 

6.87/25 Membrane Transport -2.3 



9.0 

11 

8.7 

2.8 

2.4 

12.4 

15.4 

12.1 

12.4 

3.4 

2.7 

1.7 

1.5 

1.9 

1.5 

1 

-2.5 

-2 

1 

-2.3 

2.2 
7.4 

-6.7 

-2.4 

-6.2 

-3.8 
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Table 1 Functional classification and relative fold change of K. oxytoca BAS-10 proteins differentially expressed during 
anaerobic growth in FEC (Continued) 



3 A 

34 


Leucine ABC transporter 
subunit substrate-binding 
protein LivK 


1 i\ iV 

LIVK 




5.71/40 


j.Ud/4z 


Membrane Transport 


3 n 
o.U 


5.D 


61* 


Maltose ABC transporter 

nprinlpc:mir nrnfpin 
pel I picoi 1 1 IL |JIULcMI 


MBP 


365907919 


6.88/43 


5.96/40 


Membrane Transport 


-7 


-14 


56* 


n)-n3 l^rtn^p-hi nHi nn 

periplasmic protein 


MGLB 


3763R1 367 


6.14/36 


5.42/29 


p m K r 3 n p Tr^n^nnrt 

IVICI I IkJICll IC I lul lj|JUI l 


-2.6 


-2.5 


13 


Outer membrane protein A 


OmpA 


376393752 


5.98/38 


4.98/31 


Membrane Transport 


2.7 


1 


39 


Outer membrane protein A 


OmpA 


376393752 


5.98/38 


5.5/28 


Membrane Transport 


1 


-1.6 


43 


Outer membrane protein A 


OmpA 


376393752 


5.98/38 


4.78/32 


Membrane Transport 




-1.5 


44 


Outer membrane protein A 


OmpA 


376393752 


5.98/38 


4.86/32 


Membrane Transport 




-1.6 


62 


Outer membrane protein A 


OmpA 


376393752 


5.98/38 


5.17/32 


Membrane Transport 




-2.3 




Outer membrane protein A 


umpA 


j/ojyj/jz 


j.yo/oo 


j.Uj/jz 


Membrane Transport 




1 1 


64 


Outer membrane protein A 


OmpA 


376393752 


5.98/38 


5.39/22 


Membrane Transport 




-5 


88 


Outer membrane protein A 


OmpA 


376393752 


5.98/38 


5.46/27 


Membrane Transport 




-3.7 


45 


Outer membrane protein W 


OmpW 


365909926 


6.17/23 


5.2/22 


Membrane Transport 


-4.0 


1 


18 


Hypothetical protein 
HMPREF9694_01670 
(tricarboxylic transport) 


TctC 


376400424 


8.61/35 


6.7/26 


Membrane Transport 


-6.5 


-5.0 


48* 


Hypothetical protein 
HMPREF9694_01670 
(tricarboxylic transport) 


TctC 


376400424 


8.61/35 


7.8/29 


Membrane Transport 


-4.1 


-8.0 


71* 


Hypothetical protein 
HMPREF9694_01670 
(tricarboxylic transport) 


TctC 


376400424 


8.61/35 


8.37/29 


Membrane Transport 


-10 


-6.2 


65* 


Carbamoyl phosphate 
synthase small subunit 


CPSase 


365908404 


5.79/42 


5.4/42 


Nucleotide metabolism 


-2 


-3.3 


DO 


Multifunctional nucleoside 
diphosphate kinase 


Moll/ 

i\iai\ 


vrni 171 1 

ocoy i i / 1 i 


5.55/1 5 




Nucleotide Metabolism 


— Z.D 


A ~) 

— 4.Z 


67* 


Nucleoside diphosphate 
kinase 


NdKs 


376399440 


5.55/15 


5.42/14 


Nucleotide Metabolism 


-1.8 


-4.1 


89 


Isochorismatase hydrolase 




365909686 


5.61/24 


5.36/26 


Other 


-1.6 


1 


68 


Adenosine-3'(2'),5'- 

Ulbpi lUbpi Id Lc 1 lULIcULIUdbc 


CysQ 


365908080 


5.67/28 


5.31/29 


Other 


1 


-2.1 


72 


AIIa/I h\/H rnnprnviHp 
/\ii\yi i i y u i u|-Jci UAiuc 

reductase subunit C 


AHPC 


O / UJ7JUO/ 


5 03/21 


4 95/23 


UAIUU IcUULLIUII 


—2 3 


1 


76 


Superoxide dismutase [Fe] 


FeSOD 


376402095 


5.75/21 


5.49/22 


Oxido reduction 


1 


3.1 


73 


Superoxide dismutase [Mn] 


MnSOD 


376396893 


6.23/23 


5.84/25 


Oxido reduction 


1 


-2.6 


74* 


Superoxide dismutase [Mn] 


MnSOD 


376396893 


6.23/23 


6.01/24 


Oxido reduction 


-2.8 


-6.3 


75* 


Superoxide dismutase [Mn] 


MnSOD 


376396893 


6.23/23 


5.62/24 


Oxido reduction 


-2 


-4.5 


11 


Thiol peroxidase 


TPX 


376401245 


4.67/18 


4.8/19 


Oxido reduction 


1 


-1.9 


79 


Chaperonin 


CHA10 


376396016 


5.38/10 


4.97/16 


Protein metabolism 


1 


-3.2 


82* 


Elongation factor G 


EF-G 


365907166 


5.17/77 


5.62/50 


Protein metabolism 


1.8 


1.9 


86* 


Ribosomal protein L1 


RPL1 


376396770 


9.56/25 


7.77/27 


Protein metabolism 


-2 


-2.2 


87* 


Ribosomal protein L1 


RPL1 


206564980 


9.56/25 


7.56/28 


Protein metabolism 


-1.7 


-2.1 


80 


Ribosomal protein L13 


RPL13 


365907069 


9.60/16 


8.84/16 


Protein metabolism 


1 


-3.9 


85* 


Peptidyl-prolyl cis-trans 
isomerase SurA 


SurA 


365908414 


6.42/47 


5.68/48 


Protein metabolism 


-2.1 


-2.3 
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Table 1 Functional classification and relative fold change of K. oxytoca BAS-10 proteins differentially expressed during 
anaerobic growth in FEC (Continued) 



78 Autonomous glycyl radical 
cofactor 


GrcA 


365911773 


4.82/14 


4.71/12 


Stress response 


-1.9 


1 


69 Osmotically-inducible 
protein Y 


OSMY 


376376350 


8.67/21 


5.84/22 


Stress response 


-3.0 


-3.7 


47 Universal stress protein F 


USF 


376401205 


5.46/16 


5.11/14 


Stress response 


-1.5 


1 



Fold change refers to normalized spot volumes as calculated by in silico analysis of 2D-DIGE protein maps. Positive and negative values stand for up- and down- 
regulation in anaerobic growth on FEC, respectively. Protein spots differentially abundant in anaerobic FEC with the respect to both aerobic FEC and anaerobic 
NAC are marked by an asterisk. 

functional classification was performed according to KEGG database (http://www.genome.jp/kegg/). 

b Relative fold change was measured according to the criteria reported in the experimental section. FEC: Fe(lll)-citrate containing medium. NAC: Na(l)-citrate 
containing medium. +0 2 : aerobiosis. -0 2 : anaerobiosis. 



proteolytic fragmentation. Concerning CitE and CitF, 
this is the first study suggesting PTM regulatory events. 
The differential analysis showed that anaerobic con- 
ditions promote citrate fermentation enzymes accu- 
mulation, similarly to what already reported in other 
enterobacteria [10,30,31]. In addition and more 



interestingly, during BAS-10 anaerobic growth, Fe(III) 
positively controls the accumulation of citrate fermenta- 
tive enzymes better than Na(I). In fact, in K. oxytoca 
BAS-10 citrate fermentation appears much more up- 
regulated in presence of Fe(III) than Na(I), making this 
process unique for this microorganism among 



Spot N. 



47 
69 
78 
85 
86 
87 
82 



37 

20 

19 

23 

22 

33 

41 

17 

2 

21 

28 

38 

46 

16 

15 

36 

26 

14 

81 

25 

24 

35 

32 

31 

30 

29 

27 

29 



+0 2 



Regulators 
TCA cycle 



Glycolysis/ 
gluconeogenesis 



Citrate catabolism 
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Figure 2 Functional distribution and abundance fold change (blue diamond) of differentially represented protein spots resulting from 
the proteomic comparison between anaerobic and aerobic growth on FEC. Fold change and spot number refer to Table 1. A. a. stands for 
Amino acid. Met. stands for metabolism. 
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82 
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Figure 3 Functional distribution and abundance fold change (blue diamond) of differentially represented protein spots resulting from 
the proteomic comparison between anaerobic growth on FEC and NAC. Fold change and spot number refer to Table 1. A. a. stands for 
Amino acid. Met. stands for metabolism. Cof. biosynth. stands for Cofactor biosynthesis. 



enterobacteria which usually ferment citrate in a NA(I)- 
dependent manner [1,8-10]. 

Occurrence of Fe(III) during anaerobic growth 
also modulated the abundance of many central car- 
bon metabolism enzymes. In particular, pyruvate ki- 
nase (PK), glyceraldehyde phosphate dehydrogenase 
(GAPDH) and triosephosphate isomerase (TIM) were 
up-regulated whereas TCA cycle enzymes, such as 
dihydrolipoyllysine-residue succinyltransferase (SucB), 
malate dehydrogenase (MDH) and fumarate hydratase 
(FH) were down-regulated in anaerobic FEC with re- 
spect to both aerobic FEC and anaerobic NAC (Table 1, 
Figures 2 and 3). Interestingly, fumarate reductase flavo- 
protein subunit (FrdA) was up-regulated during the an- 
aerobic growth on FEC with respect to both NAC and 
aerobic FEC (Figure 4). FrdA is part of complex II 
homolog menaquinol:fumarate oxidoreductase, which 
oxidizes menaquinol and transfers the electrons to fu- 
marate during bacterial anaerobic respiration, with 



fumarate as the terminal electron acceptor [32], thus 
counteracting TCA cycle down-regulation. Altogether 
these data suggest that to efficiently divert the carbon 
flux towards acetate and ATP production citrate fermen- 
tation enzymes are up-regulated during anaerobic 
growth on FEC, whereas TCA cycle enzymes are 
repressed (Table 1 and Figure 4). In addition, these data 
indicated that the increased catabolism of citrate 
throughout a fermentative pathway is coupled to the 
synthesis of metabolic precursors necessary for anabolic 
processes like sugar biosynthesis. In particular, the TIM 
product glycerone-P is a precursor involved in rhamnose 
biosynthesis (Figure 4). Since rhamnose is the major 
sugar of EPS [17], the observed TIM up-regulation can 
represent an interesting link between central carbon 
metabolites and EPS synthesis in Klebsiella. As a conse- 
quence of the increased anaerobic citrate fermentation, 
the production of acetic acid may also determine an in- 
crement of H + gradient across cell membrane, which 
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• Fe 3+ Up-regulation FEC Vs NAC Down-regulation FEC Vs NAC 
citrate ^ Up-regulation [-0 2 ] Vs [+0 2 ] ^ Down-regulation [-0 2 ] Vs [+0 2 ] 

Figure 4 Synoptic scheme of metabolic pathways involved in Fe(lll)-citrate catabolism. Anabolic routes, regulatory and 
membrane-associated proteins are also indicated. Reactions are reported according to KEGG [22] and EcoCy [23] databases. OM: outer 
membrane. PG: peptidoglycan. IM: inner membrane. 



positively affects the activity of the ATP synthase com- 
plex. In agreement with this view, ATP synthase subunit 
B was observed as up-regulated under anaerobic condi- 
tions in FEC medium, in the respect of both NAC and 
aerobic FEC (Table 1, Figures 2, 3 and 4). 

Central carbon and nitrogen metabolism key regulators 

The proteomic comparison showed down-regulation 
in anaerobic FEC of EIIA Glc , a component of bacter- 
ial phosphoenolpyruvate (PEP): carbohydrate phospho- 
transferase system (PTS), which was revealed as two 
protein spots differing for Mw in 2D-maps (Table 1, 
Figures 2 and 3; Additional file 1 Table SI and 
Figure S2). The differential regulation of EIIA Glc may be 



related to the differential regulation of central carbon 
metabolism enzymes. In fact, in both Gram-negative 
and Gram-positive bacteria PTS consists of several 
factors interacting with different regulatory proteins 
thus controlling glucose metabolism and many other 
cellular functions [33]. In particular, EIIA Glc is the 
central processing unit of carbon metabolism in en- 
teric bacteria since it is involved in the regulation of 
adenylate cyclase (AC) and therefore in carbon catab- 
olite repression trough the control of the catabolite 
repressor protein (CRP) (Figure 5). In addition, it 
also interacts with several non-PTS permeases 
and glycerol kinase to inhibit their activity (inducer 
exclusion) [33]. 
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Figure 5 Regulatory network controlling the expression or the activity of K. oxytoca BAS-10 differentially abundant proteins. The 

relationships are reported according to EcoCyc database [23], Bott (1997) [1], Salmon et al. (2005) [34] and Kumar et al. (201 1) [49]. Dashed lines 
refer to an indirect relationship. 



Moreover, during anaerobic growth, the down regula- 
tion of EIIA Glc may be related to the observed up- 
regulation of Arc A in FEC with respect to NAC (Table 1; 
Figures 2 and 3). ArcA, a negative response regulator of 
genes in aerobic pathways [34,35], is a central regulator 
of carbon metabolism which, from one hand, negatively 
controls the expression of pts and of TCA cycle genes 
[33,36] (Figure 5), and, from the other one, it posi- 
tively regulates the expression of genes, like citAB, pfl 
and ack, necessary for fermentative metabolic activities 
under anaerobic conditions [34,37] (Figure 5). In 
addition to its role of activator of fermentative genes 
and inhibitor of aerobic catabolism genes ArcA 
appears to regulate a wide variety of processes like 
amino acid and sugar transport and metabolism, cofac- 
tor and phospholipids biosynthesis, nucleic acid metab- 
olism [34,38]. 

Interestingly, the observed down-regulation of EIIA Glc 
in anaerobic FEC with respect to both NAC and aerobic 
FEC seemed coupled with a differential regulation of im- 
portant players of nitrogen and amino acid metabolism. 
In particular, GlnK, identified as three different spots in 
K. oxytoca BAS-10 2D-maps (Table 1, Figure 2 and 3; 
Additional file 1 Table SI and Figure S2), is part of nitrogen 
regulatory system (ntr) which facilitates the efficient assimi- 
lation of nitrogen from a variety of compounds into 



glutamine and glutamate as reported in many bac- 
teria like E. coli [39]. The activity of GlnK, which is 
encoded by an operon encoding the ammonia trans- 
porter AmtB too, is controlled by uridylylaton state. 
In nitrogen proficiency, not-uridylylated GlnK form 
is the most abundant and interacts with and inhibits 
AmtB [40]. In nitrogen limitation, uridylylated GlnK 
inhibits the phosphatase activity of the sensor mem- 
brane protein GlnL, that forms a two-component 
system with the regulatory protein GlnG. Phosphory- 
lated GlnG positively affects the expression of many ntr 
genes, including GlnK [39,41] (Figure 5). The difference in 
pi of the GlnK forms here identified gives count for protein 
uridylylation, which decreases the predicted pi value of na- 
tive protein form. The protein form with decreased Mw 
may be due to the possible removal of a N-terminal amino 
acid stretch as reported in Streptomyces coelicolor [42]. 
Interestingly, the putative uridylylated GlnK was down- 
regulated in anaerobic FEC with respect to both anaerobic 
NAC and aerobic FEC thus suggesting that nitrogen limita- 
tion was sensed by BAS-10 growing in anaerobic FEC. The 
GlnK abundance was in agreement with the down- 
regulation in anaerobic FEC of carbamoyl-phosphate syn- 
thase small subunit (CPSase) in the respect of both NAC 
and aerobic FEC. In fact, CPSase, an enzymes involved in 
pyrimidine and arginine biosynthesis, regulates nitrogen 
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disposal by catalyzing synthesis of carbamoyl phosphate 
from ammonia or glutamine and carbamate [43]. 

Carbon and nitrogen metabolism could also be asso- 
ciated to the regulation of membrane proteins participat- 
ing to amino acid up-taking, like the glutamine ABC 
transporter periplasmic-binding protein (GlnH) and glu- 
tamate and aspartate transporter subunit (DEBP), down- 
regulated in anaerobic FEC with respect to both NAC 
and aerobic FEC, or arginine 3 rd transport system 
periplasmic-binding protein (ArtJ) and leucine ABC 
transporter subunit substrate-binding protein (LivK), 
up-regulated in anaerobic FEC with respect to both 
NAC and aerobic FEC. In many bacteria, nitrogen and 
carbon metabolism are linked to cell energetic state 
through the conversion of TCA cycle intermediate 
alpha-ketoglutarate and the amino acid glutamate [44- 
46]. In K. oxytoca BAS-10 growing in anaerobic FEC, 
carbon metabolism regulators, like Arc A and CRP, and 
nitrogen assimilation factors, such as GlnG, could be 
coordinately controlled (Figure 5) in order to balance 
amino acid content and to promote biomass 
accumulation. 



Redox balance and stress response proteins 

Two isoforms of superoxide dismutase Mn-dependent 
(Mn-SOD) were down- regulated in anaerobic FEC with 
the respect to both aerobic FEC and anaerobic NAC 
(Table 1, Figure 2 and 3; Additional file 1 Table SI and 
Figure S2). Mn-SOD is devoted to destroy radicals which 
are normally produced within the cells. Interestingly and 
in agreement with what already observed in E. coli [47] 
and S. typhimurium [48], the presence of iron causes the 
accumulation of a Fe-SOD whose activity is devoted to 
contrast iron dependent free radical formation. In par- 
ticular, Fe-SOD had an abundance level not depending 
on oxygen, since their relative amount was the same 
during aerobic and anaerobic growth on FEC. Transcrip- 
tional regulation of both SOD genes is reported to be 
mediated by the Ferric uptake regulator (Fur). Fur con- 
trols the transcription of genes, like fecABCDE, involved 
in iron homeostasis and participates also in the tran- 
scriptional regulation of many other genes, such as TCA 
cycle and PTS genes [49] (Figure 5). Indeed, Fur is 
mainly a transcriptional repressor when it binds Fe(II). 
Anyway, by repressing the expression of the non-coding 
RNA Ryhb, that negatively controls mRNA translation 
into proteins, consequently Fur may also indirectly exert 
a positive gene regulatory effect as observed in E. coli 
[50]. This phenomenon could also occur in K. oxytoca 
BAS-10 since both fur and ryhb homologues were iden- 
tified in Klebsiella genus (BLAST analysis, data not 
shown). Thus, it is likely that the Fe(III) -dependent ac- 
cumulation of Fe-SOD under anaerobiosis may depend 



on Fur action to counteract metal toxicity, as observed 
in E. coli [51]. 

Furthermore, two proteins, the alkyl hydroperoxide re- 
ductase subunit C (AHPC) and the thiol peroxidase 
(TPX), contrasting reactive oxygen species (ROS), were 
differentially regulated in K. oxytoca BAS-10. In particu- 
lar, AHPC was down-regulated in anaerobic FEC 
with respect to aerobic FEC, thus revealing an oxygen 
dependent control, while TPX was down-regulated in 
anaerobic FEC with respect to anaerobic NAC, thus re- 
vealing a Fe(III) -dependent regulation. 

Concerning stress associate proteins, the hyperosmoti- 
cally inducible periplasmic (OSMY) was down- regulated 
in anaerobic FEC with respect to both aerobic FEC and 
anaerobic NAC. This protein was found associated with 
osmotic stress response during E. coli [52]. Other stress 
associated proteins, like DNA protection during starva- 
tion protein (DPS) and the universal stress protein F 
(USF), were down-regulated in anaerobic FEC with the 
respect to aerobic FEC showing an oxygen dependent 
regulation. In DPS, a ferritin-like protein, binds to the 
chromosome and protects DNA from oxidative damage 
by sequestering intracellular Fe(III) and storing it in the 
form of Fe(III)-oxyhydroxide mineral [53]. 

Proteome data revealed that outer membrane proteins 
are differentially regulated under anaerobic condition in 
presence of FEC or NAC. In particular, three isoforms of 
major outer membrane murein-associated lipoprotein 
(Lpp), differing in pi value, were up-regulated in anaer- 
obic FEC. In E. coli, Lpp is synthesized as a precursor 
protein (prolipoprotein) and processed by a signal pep- 
tidase after modification with the addition of di-acyl- 
glycerol and a fatty acid chain [54]. No modification 
altering pi has been reported up to day being this study 
the first one suggesting the possibility of such similar 
PTMs in Lpp. This finding was coupled with the down- 
regulation of outer membrane OmpA. Indeed, OmpA 
was revealed as four down- regulated full size isoforms 
and one and three protein fragments up- and down- 
regulated in anaerobic FEC, respectively (Table 1). In E. 
coli OmpA acts as a low permeability porin and it is 
present in E. coli protein 2D -maps as different protein 
spots suggesting PTM control [55]. Interestingly, both 
Lpp and OmpA are reported to interact with PAL that is 
part of Tol-PAL complex, a membrane-spanning multi- 
protein system that has been reported playing several 
functions in Gram-negative bacteria, including transport 
regulation, cell envelope integrity and pathogenicity [56]. 
In addition, Lpp has been reported to interact with TonB 
protein, which serves to couple the cytoplasmic mem- 
brane proton motive force to the active transport of 
iron-siderophore complexes, including the Fee system 
[57]. The differential regulation of Lpp and OmpA sug- 
gests differential roles for these two proteins in 
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controlling iron omeostasis, membrane transport and/or 
preventing possible cellular damage. 

In conclusion, this proteomic analysis revealed an 
iron-dependent regulation for several stress-related pro- 
teins; altogether, these observations suggest the occur- 
rence of different concomitant mechanisms to contrast 
and/or prevent BAS-10 cell damage during aerobic or 
anaerobic growth. 

Conclusion 

At the best of our knowledge, for the first time this 
study describes at the proteome level biochemical 
mechanisms allowing the physiological adaptation of the 
enter obacterium K. oxytoca strain BAS-10 [12,16-18] to 
sustain anaerobic growth on Fe (III) -citrate as sole car- 
bon and energy source. K. oxytoca BAS-10, isolated from 
pyrite mine drainages of Colline Metallifere (Tuscany, 
Italy) and sharing the same ecological niche of specia- 
lized bacteria like Acidithiobacillus ferrooxidans [58], has 
the peculiarity of thriving on high concentrations of Fe 
(III) -citrate. This character distinguishes it from other 
clinically-isolated enterobacteria, belonging to K. oxy- 
toca, K. pneumonia, S. typhimurium and Neisseria 
gonorrhoeae species, which are able to uptake and fer- 
ment anaerobically citrate as sole carbon and energy 
source in a Na(I) -dependent manner only [1,10,12,59]. 
When anaerobically grown on FEC or NAC, K. oxytoca 
BAS-10 biomass production yields and citrate consump- 
tions were very similar [16]. Under this condition, about 
half of the initial Fe(III) was reduced to Fe(II) [12,16], 
the highest levels observed so far with respect to other 
fermentation processes [60,61]. 

Differential proteome analyses, carried out using an- 
aerobic FEC as pivotal condition, revealed two sets of 
proteins whose abundance is 0 2 - and Fe(III) -dependent, 
respectively. The cross checking of these proteome data 
revealed proteins whose abundance is specifically asso- 
ciated to anaerobic growth on Fe(III). Interestingly, these 
data portrayed a coordinated regulation of citrate cata- 
bolic enzymes during anaerobic growth on FEC. Citrate 
may enter the cell throughout several kind of permeases, 
including Cit transporters [1,10], TctABC [3] and 
FecABCDE systems [4-7]. Then, it is catabolised via the 
TCA cycle or the citrate fermentation pathway under 
aerobic or anaerobic growth, respectively. Under anaer- 
obic conditions, the first pathway is repressed, while the 
second one is favoured. In addition, a fascinating specu- 
lation is that one BAS-10 citrate fermentation enzyme, 
the oxalacetate decarboxylase, can be adapted to act as 
both Fe(II)- and Na(I) -depending pump to excrete Fe(II) 
and Na(I). This hypothesis has to be verified by func- 
tional studies on BAS-10 oadGAB genes. Furthermore, 
PK, GAPDH and TIM abundance levels suggested an 
increased synthesis of metabolic precursors necessary 



for anabolic processes, such as those yielding rhamnose, 
linking together central carbon metabolism and EPS bio- 
synthesis in Klebsiella. Differential regulation of meta- 
bolic enzymes was associated with variable cellular levels 
of carbon metabolism key regulators ArcA and PTS 
EIIA Glc , and important players of nitrogen and amino 
acid metabolism, like GlnK. In this context, ArcA plays 
a central role in regulating anaerobic carbon metabolism 
by controlling the expression of citrate fermentation and 
TCA cycle enzymes directly or indirectly - i.e. through- 
out PTS gene regulation - and nitrogen metabolism (Fig- 
ure 5). Furthermore, during anaerobic FEC growth, 
ArcA-dependent gene regulation may be combined to 
iron effect on gene expression trough the activity of Fur 
(Figure 5). Proteome results also suggest that metabolic 
pathway modulation seems coupled with regulation of 
redox balance, stress response and outer membrane pro- 
tein content. The occurrence of various proteins involved 
in these processes is finely balanced to face iron and pH 
toxicity, thus allowing K. oxytoca BAS-10 to survive in 
extreme habitats. Altogether, these data well correlate 
with previous observations on other microorganisms able 
surviving in highly polluted environments [62,63] or ex- 
tremophile organisms living under extreme conditions 
[64,65], both generally activating multiple intracellular 
processes devoted to contrast environmental toxic effects 
and optimize energy production. 

According to the results presented here, K. oxytoca 
BAS-10 faces high iron concentration by various mo- 
lecular adaptation mechanisms, which are summarized 
in Figures 4 and 5. The coordinated regulation of meta- 
bolic pathways and cellular networks allows K. oxytoca 
BAS-10 to efficiently sustain bacterial growth by adapt- 
ing metabolic and biochemical processes in order to face 
iron toxicity and to ensure optimization of energy pro- 
duction. In this context, EPS production is intrinsically 
related to bacterial biomass production and metabolic 
and energy supply for anabolic reactions. Since K. oxy- 
toca BAS-10 was proven to have a high potential for bio- 
technological applications [18-20], the data presented 
here could be used for redesigning fermentation strat- 
egies that are difficult to be intuitively identified and 
approaching novel genetic targets to be engineered for a 
rational design of high-yielding EPS producer strains. 

Methods 

Bacterial cultivation and metal binding exopolysaccharide 
(EPS) quantification 

For K. oxytoca BAS-10 cultures, two defined media 
were prepared containing the same mineral composition 
(2.5 g NaHC0 3 , 1.5 g NH 4 C1, 1.5 g MgS0 4 .7H 2 0, 0.6 g 
NaH 2 P0 4 , 0.1 g KC1, buffered at pH 7.6 with NaOH) 
and 50 mM Na-citrate (14.7 g.l" 1 ), named NAC medium, 
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or 50 mM ferric citrate (13.15 g.L 1 ), named FEC 
medium, as the sole carbon and energy source. BAS-10 
was retrieved from cryovials kept at -80°C in 25% gly- 
cerol in Nutrient Broth (NB) (Difco). One ml aliquots of 
dense (abs 60 onm = 1.4) overnight NB pre-culture of BAS- 
10 were used to inoculate (1:100 v:v) FEC or NAC 
medium cultivations. Aerobic cultivations were per- 
formed in flasks, aerated by a magnetic stir bar (250 r.p. 
m.). Anaerobic cultivations were performed in pirex bot- 
tles, previously fluxed with N 2 and kept under anaero- 
biosis by sealed cap. Aerobic and anaerobic cultivations 
were performed in parallel quadruplicates, at 30°C. Bio- 
mass samples for proteomic analysis were harvested by 
centrifugation at late exponential growth stages, i.e. after 
36 and 72 h of incubation in aerobic and anaerobic culti- 
vations, respectively. For EPS extraction and quantifica- 
tion, the procedures described by Baldi et al. (2009) [12] 
were followed. 

Phylogenetic analysis 

One microliter of diluted (1:10 in autoclaved distilled water) 
BAS-10 bacterial suspension was used for PCR amplifica- 
tion of 16S rRNA gene by using fDl and rDl universal pri- 
mers [66] and 1.5 units of recombinant Taq DNA 
Polymerase (Invitrogen, Life Technologies), according to 
manufacturer s specifications. To perform bacterial cell lysis 
and DNA denaturation, a treatment of 5 min, at 95 °C, was 
performed. Then, amplification steps were carried out for 
40 cycles consisting of 1 min at 95°C, 1 min at 55°C and 2 
min at 72 °C, with a final extension of 10 min at 72 °C. 
PCR fragments were purified by using PCR clean-up kit 
(NucleoSpin, Macherey-Nagel) according to manufacturer s 
specifications and sequenced with the same primers at 
BMR Genomics (University of Padova, Italy). 

BAS-10 16S rDNA sequence was compared using 
BLAST probing [67] of DNA sequences from the 
NCBI database (http://blast.ncbi.nlm.nih.gov/Blast.cgi? 
CMD=Web&PAGE_TYPE=BlastHome) with default pa- 
rameters, selecting either only cultivable reference 
strains (refseq_rna) or whole database strains (rn/nt). 
ClustalX program [68] was used to align the homolo- 
gous 16S rDNA gene sequences obtained from the data- 
base, choosing the slow/accurate method for pairwise 
alignment. All positions containing alignment gaps were 
eliminated in pairwise sequence comparisons by select- 
ing the related option. The 16S rDNA sequence of Strep- 
tomyces tendae was used as outgroup. Alienated 
sequences were then used to generate rooted trees 
throughout Unweighted Pair Group Method with Arith- 
metic Mean (UPGMA) algorithm. 

Protein extraction 

Total proteins were extracted from frozen biomass sam- 
ples of quadruplicated cultivations per each condition by 



using an experimental procedure previously described 
[69,70]. Briefly, pellets were cleaned three times with a 
washing solution (10 mM Tris-HCl pH 7.5, 5 mM 
EDTA, 1 mM DTT, 0.5 mM PMSF, 4 mg/ml leucopep- 
tin, 0.7 mg/ml pepstatin, 5 mg/ml benzamidin). Follow- 
ing re-suspension in washing solution containing 0.3% 
SDS, bacterial cells were disrupted by sonication on ice 
(output control 4, 4 x 15 s, Vibra Cell, USA). Samples 
were then boiled (5 min) and rapidly cooled down on ice 
(15 min). DNase (100 jig/ml) and RNase (50 jig/ml) 
were added in ice for 20 min. Cell debris and non- 
broken cells were separated by centrifugation at 15,000 x g, 
for 15 min, at 4°C. Protein solutions were dialyzed against 
distilled water for 3 h, at 4°C, precipitated by using 3 vol 
of acetone, at -20°C, overnight, and re-suspended in 
30 mM Tris pH 8.5, 7 M urea, 2 M thiourea, 4% w/v 
CHAPS. 

2D-DIGE analysis 

Proteomes of K. oxytoca BAS-10 grown in FEC medium 
under anaerobiosis [-0 2 ] and aerobiosis [+0 2 ] or in 
NAC under [-0 2 ] conditions were compared, using FEC 
medium under [-0 2 ] as pivotal condition. To this aim, 
protein samples were labelled for 2D-DIGE analysis 
using CyDye™ DIGE minimal labelling kit (GE Health- 
care, Sweden), as previously described [71]. Briefly, pro- 
tein samples (40 \ig) were labelled with 320 pmol of 
CyDye on ice in the dark for 30 min. Labelling reaction 
was stopped by addition of 0.8 \A of 10 mM lysine and 
incubation was continued on ice for 15 min., in the dark. 
Two samples out of a total of four per each condition 
were labelled using Cy3 dye and the other two using 
Cy5 dye to account for florescence bias. In addition, a 
standard protein pool was generated by combining an 
equal amount of protein extracts from which six 40 \ig 
protein aliquots were then minimally labelled with 
320 pmol Cy2 fluorescent dye (CyDye™ DIGE, GE 
Healthcare), according to manufacturers instructions. 
Thus, a total of six 2D-DIGE gels were performed con- 
taining a mix of Cy2-labeled pooled protein standard 
and Cy3- and Cy5-labeled protein samples from cells 
incubated in FEC medium under [-0 2 ] and [+0 2 ] and 
in NAC under [-0 2 ] conditions, analysed in couples, i.e. 
FEC [-0 2 ] vs FEC [+0 2 ], FEC [-0 2 ] vs NAC [-0 2 ] and 
FEC [+0 2 ]vs NAC [-0 2 ]. For IEF, DeStreak rehydration 
solution (GE Healthcare), containing 0.68% v/v IPG buf- 
fer (GE Healthcare) and 10 mM DTT (Sigma), was 
added to each mix up to a 340 \A final volume. IEF was 
performed as previously described [57] using 4-7 pH 
range 18 cm-IPG strips (GE Healthcare) in an Ettan 
IPGphor III apparatus (GE Healthcare). After IEF, IPG 
strips were incubated with an equilibration buffer (6 M 
urea, 30% v/v glycerol, 2% w/v SDS, 0.05 M Tris-HCl, 
pH 6.8) containing 2% w/v DTE for 10 min; thiol groups 
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were then blocked by further incubation with equilibra- 
tion buffer containing 2.5% w/v iodoacetamide. Focused 
proteins were then separated by using 12% SDS-PAGE, 
at 10°C, in a Hettan Dalt six (GE Healthcare), with a 
maximum setting of 40 \iA per gel and 110 V. 

Protein visualization and image analysis 

The 2D-gels were scanned with a DIGE imager (GE 
Healthcare) to detect cyanin-labeled proteins, according 
to manufacturers instructions. Differential gel analysis 
was performed automatically by using the Image Master 
2D Platinum 7.0 DIGE software (GE Healthcare), 
according to the manufacturers instructions. Protein 
spots were automatically detected and then matched. In- 
dividual spot abundance was automatically calculated 
from quadruplicated 2D-gels as mean spot volume, i.e. 
integration of optical density over spot area, and normal- 
ized to the Cy2-labeled internal pooled protein standard. 
Protein spots showing more than 1.5 fold change in spot 
volume (increased for up-regulation or decreased for 
down-regulation), with a statistically significant ANOVA 
value (P < 0.05), were considered differentially repre- 
sented and further identified by MS analysis. 

Protein identification 

Protein spots were excised from the 2D-gels, alkylated, 
digested with trypsin and identified as previously 
reported [56]. Peptide mixtures were desalted by uZip- 
TipC18 (Millipore, MA) using 50% v/v acetonitrile/5% 
v/v formic acid as eluent before MALDI-TOF-MS and 
nLC-ESI-LIT-MS/MS analysis. 

In the case of MALDI-TOF-MS experiments, peptide 
mixtures were loaded on the MALDI target, by using 
the dried droplet technique, a-cyano-4-hydroxycinnamic 
acid as matrix, and analyzed using a Voyager-DE PRO 
mass spectrometer (Applied Biosystems, USA) operating 
in positive ion reflectron, with an acceleration voltage of 
20 kV, a N 2 laser (337 nm) and a laser repetition rate of 
4 Hz [72]. Final mass spectra, measured over a mass 
range of 700-4000 Da and obtained by averaging 400- 
800 laser shots, were elaborated using the DataExplorer 
5.1 software (Applied Biosystems) and manually 
inspected to get the corresponding peak lists. Internal 
mass calibration was performed with peptides deriving 
from trypsin autoproteolysis. Tryptic digests were even- 
tually analyzed by nLC-ESI-LIT-MS/MS using a LTQ XL 
mass spectrometer (Thermo, San Jose, CA) equipped 
with a Proxeon nanospray source connected to an Easy- 
nanoLC (Proxeon, Odense, Denmark) [73,74]. Pep- 
tide mixtures were separated on an Easy CI 8 column 
(10-0.075 mm, 3 \im) (Proxeon). Mobile phases were 
0.1% v/v aqueous formic acid (solvent A) and 0.1% v/v 
formic acid in acetonitrile (solvent B), running at total 
flow rate of 300 nL/min. Linear gradient was initiated 



20 min after sample loading; solvent B ramped from 5% 
to 35% over 15 min, from 35% to 95% over 2 min. Spec- 
tra were acquired in the range m/z 400-1800. Acquisi- 
tion was controlled by a data-dependent product ion 
scanning procedure over the 3 most abundant ions, en- 
abling dynamic exclusion (repeat count 2 and exclusion 
duration 60 s); the mass isolation window and collision 
energy were set to m/z 3 and 35%, respectively. 

MASCOT search engine version 2.2.06 (Matrix Sci- 
ence, UK) was used to identify protein spots from an 
updated NCBI non-redundant database by using 
MALDI-TOF-MS and nLC-ESI-LIT-MS/MS data and by 
selecting trypsin as proteolytic enzyme, a missed clea- 
vages maximum value of 2, Cys carbamidomethylation 
and Met oxidation as fixed and variable modification, re- 
spectively. In the first case, a mass tolerance value of 
40-80 ppm was selected; in the second case, a mass tol- 
erance value of 2 Da for precursor ion and 0.8 Da for 
MS/MS fragments was chosen. MALDI-TOF-PMF can- 
didates with a MASCOT score > 83 or nLC-ESI-LIT- 
MS/MS candidates with more than 2 assigned peptides 
with an individual MASCOT score > 25, both corre- 
sponding to P < 0.05 for a significant identification, were 
further evaluated by the comparison with their calcu- 
lated mass and pi values, using the experimental data 
obtained from 2D-DIGE. 

Additional file 



Additional file 1: Table SI: containing mass spectrometry 
parameters of K. oxytoca BAS-10 differentially abundant protein 
identification; Figure SI: showing phylogenetic tree generated by 
using 16S rDNA sequence of K. oxytoca BAS-10 and the first twenty 
hits from a BLAST analysis performed by selecting whole database 
strains; Figure S2 showing 2D-protein maps, chosen as examples of 
anaerobic FEC, anaerobic NAC and aerobic FEC condition, 
respectively. Labels indicating differentially abundant protein spots and 
referring to Table 1 and Table 1 are also reported in Figure S2; Figure 
S3: containing distribution into functional classes of the identified 
protein spots, according to KEGG metabolic database (http://www. 
genome.jp/kegg/), from the comparison between anaerobic and aerobic 
growth on FEC from the comparison between anaerobic growth on FEC 
and NAC, respectively; 1447 nt gene sequence generated from BAS-10 
16S rDNA. 
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